Gravitational waves and neutrino emission from the merger of binary neutron stars 
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Numerical simulations for the merger of binary neutron stars are performed in full general rela- 
tivity incorporating a finite-temperature (Shen's) equation of state (EOS) and neutrino cooling for 
the first time. It is found that for this stiff EOS, a hypermassive neutron star (HMNS) with a long 
lifetime 10 ms) is the outcome for the total mass ^ S.OMq. It is shown that the typical total 
neutrino luminosity of the HMNS is 3-8 x 10^'^ erg/s and the effective amplitude of gravitational 
waves from the HMNS is 4-6 x 10"^^ at / = 2.1-2.5 kHz for a source distance of 100 Mpc. We also 
present the neutrino luminosity curve when a black hole is formed for the first time. 
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Introduction: Coalescence of binary neutron stars 
(BNS) is one of the most promising sources for next- 
generation kilo-meter-size gravitational-wave (GW) de- 
tectors and also a possible candidate for the progen- 
itor of short-hard gamma-ray bursts (SGRB) Moti- 
vated by these facts, numerical simulations have been ex- 
tensively performed for the merger of BNS in the frame- 
work of full general relativity in the past decade since the 
first success in 2000 [3] (see, e.g., jif for a review). 

BNS evolve due to gravitational radiation reaction and 
eventually merge. Before the merger sets in, each neutron 
star is cold (i.e., thermal energy of constituent nucleons 
is much smaller than the Fermi energy) , because thermal 
energy inside the neutron stars is significantly reduced 
by neutrino and photon coolings due to the long-term 
evolution (typically ^ 10® yrs) until the merger 5]. By 
contrast, after the merger, shocks are generated by hy- 
drodynamic interactions. In particular, when a hyper- 
massive neutron star (HMNS) is formed, spiral arms are 
developed in its envelope and continuous heating occurs 
due to the collision between the HMNS and spiral arms 
(e.g., 043)- t^i^ process the maximum temperature 
increases to ~ 30-50 MeV, and hence, copious neutrinos 
are emitted [§-[lH. To accurately study the evolution of 
the hot HMNS with a physical modeling, we have to in- 
corporate microphysical processes such as neutrino emis- 
sion and equation of state (EOS) based on a theory for 
the high-density and high-temperature nuclear matter. 
However, such simulations have not been done yet in full 
general relativity (but see for a work in an approxi- 
mate general relativistic gravity with finite-temperature 
EOSs). Incorporation of microphysical processes is in 
particular important for exploring the merger hypothesis 
of SGRB because it may be driven through pair annihi- 
lation of neutrino-antineutrino pairs 0- 

In this letter, we present the first results of numerical- 
relativity simulation for the BNS merger performed in- 
corporating both a finite-temperature (Shen's) EOS [13] 
and neutrino cooling [l^ . In the following, we report the 
possible outcome formed after the merger, criterion for 
the formation of HMNS and black hole (BH), thermal 
properties of the HMNS, and neutrino luminosity from 
the HMNS and in the BH formation. 



Setting of numerical simulations: Numerical simula- 
tions in full general relativity are performed using the 
following formulation and numerical schemes. Einstein's 
evolution equations are solved in the so-called BSSN- 
puncture formulation [l5l |. We employ the original ver- 
sion of the geometric variables for the BSSN formula- 
tion together with an improved definition of the confor- 
mal factor. As in [?!, the dynamical gauge condition for 
the lapse function and shift vector is chosen; a fourth- 
order-accurate finite differencing in space and a fourth- 
order Runge-Kutta time integration are used; a conserva- 
tive shock capturing scheme with third-order accuracy in 
space and fourth-order accuracy in time is employed for 
solving hydrodynamic equations. In addition to the ordi- 
nary hydrodynamic equations, we solve evolution equa- 
tions for the neutrino (Y^), electron {Ye), and total lep- 
ton (Yi) fractions per baryon, and, taking into account 
weak interaction processes PJ]. Shen's EOS, tabulated 
in terms of the rest-mass density (p), temperature (T), 
and Ye or Y/, is employed. In addition, we incorporate 
neutrino cooling, employing a general relativistic leakage 
scheme p^ . In our leakage scheme, electron neutrinos 
(ve), electron antineutrinos (j>e), and other types (yu/r) 
of neutrinos (vx) are taken into account. 

Numerical simulations are performed preparing a non- 
uniform grid as in ^7]. The inner domain is composed of 
a finer uniform grid and the outer domain of a coarser 
nonuniform grid. The grid resolution in the inner zone is 
chosen so that the major diameter of each neutron star in 
the inspiral orbit is covered by 60 and 80 grid points for 
low- and high-resolution runs, respectively: We always 
performed simulations for both grid resolutions to con- 
firm that convergence is approximately achieved. Outer 
boundaries are located in a local wave zone (at « 560- 
600 km along each coordinate axis which is longer than 
gravitational wavelength in the inspiral phase). During 
the simulations, we check the conservation of baryon rest- 
mass, total gravitational mass (Arnowitt-Deser-Misner 
mass plus radiated energy of GWs), and total angular 
momentum (including that radiated by GWs), and find 
that the errors are within 0.5%, 1%, and 3%, respectively, 
for the high-resolution runs within the duration « 30 ms. 

Shen's EOS, derived from a relativistic mean-field the- 
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FIG. 1: Maximum rest-mass density, maximum matter tem- 
perature, and total neutrino luminosity as functions of time 
for three models. The dashed vertical line shows the time at 
which a BH is formed for model H. 

ory [isj, is a stiff one among other EOSs, giving the max- 
imum mass of zero-temperature spherical neutron stars 
■^max ~ 2.2Mq. The latest discovery of a high-mass neu- 
tron star with mass 1.97±0.04Mq [16| suggests that stiff 
EOSs are favored, and Shen's EOS satisfies this require- 
ment. There are two possible fates of BNS [i^: If its total 
mass M is larger than a critical mass M^, a BH will be 
formed soon after the onset of the merger, while a differ- 
entially rotating HMNS wiU be formed for M < M^. The 
value of Mc depends strongly on the EOS. Because Shen's 
EOS is quite stiff, Mc is much larger than the typical to- 
tal mass of BNS, ~ 2.1 Mq [H, as shown below. Thus, 
with this EOS, a HMNS is the frequent outcome, as in 
the cases of stiff EOSs with which M^ax > 2Mq [3- 

This paper focuses on the merger of equal-mass BNS 
with three masses for each neutron star: Mns = 1.35, 1.5, 
and 1.6Mq (Mns is the gravitational mass of a neutron 
star in isolation). We refer to each model as models L 
(light), M (middle), and H (heavy). We perform the 
simulation with the initial condition of about 3-4 orbits 
before the onset of the merger until the system relaxes to 
a quasistationary state. Quasiequilibrium states of BNS 
are prepared as the initial conditions, as in 01- 

Numerical results: Figure [1] plots the evolution of 
the maximum rest-mass density, Pmax, maximum mat- 
ter temperature, Tmax, and total neutrino luminosity as 
functions oi t — tmcrgo where tmerge is the onset time of 
the merger. For t < tmorgc, Pmax is approximately con- 
stant, while t ^ tmcrgo, a HMNS is formed and subse- 
quently contracts by emission of GWs, which carry en- 
ergy and angular momentum from the HMNS; Pmax in- 
creases in the gravitational radiation timescale. However, 
at i — imcrgc ^ 20 ms for models L and M, the degree 
of its nonaxial symmetry becomes low enough that the 
emissivity of GWs is significantly reduced. Because no 
dissipation process except for neutrino cooling is present, 
the HMNS will be alive at least for the cooling timescale 
before collapsing to a BH (see below). For model H, the 
HMNS eventually collapses to a BH after the gradual 
contraction due to the GW emission and a massive disk 



of « 0.1 Mq is formed around the BH. 

The evolution of Tmax plotted in Fig. [ifb) shows that 
HMNS just after the formation are hot with T^ax ~ SO- 
TO MeV. Such high temperature is achieved due to the 
liberation of kinetic energy of the orbital motion at the 
collision of two neutron stars. Subsequently, Tmax de- 
creases due to the neutrino cooling, with the maximum 
luminosity 3-10 x lO^'^ erg/s (see Fig. [Tfc)), but relaxes 
to a high value with 25-50 MeV. Around the HMNS, 
spiral arms are formed and shock heating continuously 
occurs when the spiral arms hit the HMNS (cf. Fig. [2]). 
Due to this process and because of a long neutrino cool- 
ing timescale (see below), the temperature (and thermal 
energy) does not significantly decrease in 10-100 ms. 

Figure[2]plots the color maps of rest-mass density, mat- 
ter temperature, and total neutrino luminosity for model 
M at Emerge = 15 ms, at which it relaxes to a semi-final 
quasisteady state. This shows that the HMNS is weakly 
spheroidal and the temperature is high in its outer re- 
gion. The neutrino luminosity is also high in its outer 
region, in particular, near the polar surface. With the 
fact that the rest-mass density is relatively small near 
the rotation axis above the polar surface, this is a fa- 
vorable feature for the merger hypothesis of SGRB; pair 
annihilation of neutrinos and antineutrinos could supply 
a large amount of thermal energy which may drive a fire 
ball along the rotation axis. The pair annihilation effi- 
ciency has been approximately estimated in the previous 
works [gl-fllj. These show that the efficiency for j/ej>e an- 
nihilation is 10~^(Ljy^/ 10^'^ erg/s). If this result holds 
in our work, the pair annihilation luminosity would be 
~ 10^^ erg/s. 

The reasons that HMNS are formed are; (i) it is rapidly 
rotating with the period ^ 1 ms, and hence, the centrifu- 
gal force increases the possible mass that can be sus- 
tained; (ii) because it is hot, thermal energy enhances 
the pressure. We find that the rotational velocity with 
the period 1 ms does not play a substantial role. Ex- 
ploring in detail Shen's EOS for high density tells us that 
the effect of the thermal energy is significant and can in- 
crease Minax by ~ 20-30% for a high-temperature state 
with T ^ 20 MeV. Therefore, the HMNS will alive before 
collapsing to a BH for a long cooling time Eth/ ^ 2-3 s 
where Eth is total thermal energy of the HMNS. At the 
time when the HMNS collapses to a BH, it will be close 
to a spherical configuration with low temperature due to 
longterm GW and neutrino emissions. Thus, observable 
signals from the late-time collapse will be weak. 

Figure [3] plots neutrino luminosities as functions of 
time for three fiavors (I'e, i^e, and sum of Vx)- It is found 
that electron antineutrinos are dominantly emitted for 
any model. The reason for this is as follows: The HMNS 
has a high temperature, and hence, electron-positron 
pairs are efficiently produced from thermal photons, in 
particular in its envelope. Neutrons efficiently capture 
the positrons to emit antineutrinos whereas electrons are 
not captured by protons as frequently as positrons be- 
cause the proton fraction is much smaller. Such hierarchy 
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FIG. 2: Color maps of rest-mass density (with velocity fields), temperature, and total neutrino luminosity at t ~ f5 ms after 
the merger for model M. The upper and lower panels show the configuration in the x-y and x-z planes, respectively. 




FIG. 3: Neutrino luminosities for three flavors for three mod- 
els. The inset of the bottom panel focuses on the luminosities 
in the BH formation. The meaning of the dashed line is the 
same as in Fig. [T] 



in the neutrino luminosities was reported also in [9j. 

Soon after the BH formation for model H, /x/r neutrino 
luminosity steeply decreases because high temperature 
regions are swallowed into the BH, while luminosities of 
electron neutrinos and antineutrinos decrease only grad- 
ually because these neutrinos are emitted via charged- 
current processes from the massive accretion disk. We 
here note that magnetic fields, which are not taken into 
account in the present simulations, could be amplified 
significantly in the accretion disk Q and may play a role 
in the late evolution of the BH-disk system. 

The antineutrino luminosity for the long-lived HMNS 
is Lp ~ 1.5-3 X 10^^ erg/s with small time variabil- 
ity. It is by a factor of ~ 1-5 larger than that from 



protoneutron stars formed after supernovae [19| . Aver- 
aged neutrino energy is ep ~ 20-30 MeV. The sensitiv- 
ity of water-Cherenkov neutrino detectors such as Super- 
Kamiokande and future Hyper-Kamiokande (HK) have a 
good sensitivity for such high-energy neutrinos in partic- 
ular for electron antineutrinos |20j. The detection num- 
ber for electron antineutrinos is approximately estimated 
by aATL^/{ATTD^e^) where a is the total cross section 
of the detector against target neutrinos, AT is the life- 
time of the HMNS, and D is the distance to the HMNS. 
For a one-Mton detector such as HK, the expected de- 
tection number is ^ 10 for D <^ 5 Mpc with AT ~ 2-3 s, 
based on an analysis of ;20]. Thus, if the BNS merger 
fortunately happens within Z? ~ 5 Mpc, neutrinos from 
the HMNS may be detected and its formation may be 
confirmed. Note that GWs from the HMNS will be si- 
multaneously detected for such a close event (see below) , 
reinforcing the confirmation of the HMNS formation. 

Figure|4lja) plots gravitational waveforms as a function 
of retarded time t^et ^ t - D ~ 2Mlog{D/M) for three 
models where M — 2AfNS- Here, h+ and denote the 
plus and cross modes of GWs extracted from the metric 
in the local wave zone. The waveforms are composed of 
the so-called chirp waveform, which is emitted when the 
BNS is in an inspiral motion (for tict ^ tmcrgo), and the 
merger waveform (for t^ct ^ ^morgc), on which we here 
focus. For the HMNS formation, the merger waveforms 
are composed of quasiperiodic waves for which h < 10^^^ 
for D — 100 Mpc and the peak frequencies are in a nar- 
row range /peak — 2.1-2.5 kHz depending weakly on AI. 
They agree with that in the approximate general rela- 
tivistic study [l^]. Note that /peak depends on adopted 
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FIG. 4: (a) GWs observed along the axis perpendicular to the orbital plane for the hypothetical distance to the source 
D — 100 Mpc. (b) The effective amplitude of GWs as a function of frequency for D — 100 Mpc. The noise amplitudes of 
Advanced Laser Interferometer Gravitational wave Observatories (adv. LIGO), broadband configuration of Advanced LIGO 
(bro. LIGO), and Large-scale Cryogenic Gravitational wave Telescope (LCGT) are shown together. 



EOS [T2|, and we will describe the dependence of GWs 
on EOS elsewhere. The accumulated effective amplitude, 
/icff = 0.4/i(/AT)^/^, is much larger where the factor 0.4 
comes from the averages of angular direction of the source 
and rotational axis of the HMNS. Figure |4l[b) shows the 
effective amplitude defined by 0.4ft(/)/ ^ 4-6 x 10^^^ 
for D ~ 100 Mpc, where h{f) is the absolute value of the 
Fourier transformation of /i+ + ihx- This suggests that 
for a specially-designed version of advanced GW detec- 
tors such as broadband LIGO, which has a good sensi- 
tivity for a high-frequency band, GWs from the HMNS 
oscillations may be detected with S /N= 5 if £> ^ 20 Mpc 
or the source is located in an optimistic direction. 

Summary: We have reported the first results of the 
numerical-relativity simulation performed incorporating 
both a finite-temperature (Shen's) EOS and neutrino 
cooling effect. We showed that for such a stiff EOS, 
HMNS is the canonical outcome and BH is not promptly 
formed after the onset of the merger as long as the total 
mass of the system is smaller than 3.2Mq. The primary 



reason is that thermal pressure plays an important role 
for sustaining the HMNS. We further showed that the 
lifetime of the formed HMNS with mass ^ 3Mq is much 
longer than its dynamical timescale, 3> 10 ms, and will 
be determined by the timescale of neutrino cooling. Neu- 
trino luminosity of the HMNS was shown to be high as 
~ 3- 10 X 10^^ erg/s. The effective amplitude of GWs is 
4-6 X 10~^^ at /peak = 2.1-2.5 kHz for a source distance 
of 100 Mpc. If the BNS merger happens at a relatively 
short source distance or is located in an optimistic direc- 
tion, such GWs may be detected and HMNS formation 
will be confirmed. 
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